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ABSTRACT 

Using Spitzer IRAC observations from the SAGE-SMC Legacy program and archived Spitzer IRAC 
data, we investigate dust production in 47Tuc, a nearby massive Galactic globular cluster. A previous 
study detected infrared excess, indicative of circumstellar dust, in a large population of stars in 47Tuc, 
spanning the entire Red Giant Branch (RGB). We show that those results suffered from effects caused 
by stellar blending and imaging artifacts and that it is likely that no stars below «1 mag from the tip 
of the RGB are producing dust. The only stars that appear to harbor dust are variable stars, which 
are also the coolest and most luminous stars in the cluster. 

Subject headings: circumstellar matter — globular clusters: individual (47Tuc) — infrared: stars — 
stars: AGB and post-AGB — stars: mass-loss — stars: winds, outflows 



1. INTRODUCTION 

Stellar mass loss is critical to the late stages of stel- 
lar evolution. Mass loss in Red Giant Branch (RGB) 
stars determi nes cluster Horizontal Bran ch (HB) mor- 
phology (e.g., lRoolll97l iCatelanl 12000) and the pul- 
sational properties of RR Lyr ae variables (e.g.. iChristvl 
H966tlCaloi fc D'Antonal 12008ft . Strong mass loss on the 
Asymptotic Giant Branch (AGB) can exceed the nuclear 
consumption rate, de termining the timescale of subse- 
quent evolution (e.g.. Ivan Loon et al.l 119991 : iBover et al.1 
I2009H ). AGB winds are enriched due to a convective pro- 
cess that brings processed material to the stella r surface 
(the third dredge-up; e.g.. iRenzini fc Volilll98lD . making 
AGB stars a major source of galaxy enrichment (|Gehrzl 
[1989] ). In the low -mass stars of old globular clusters 
(GCs; M < 0.8 M R ), more mass is los t on the RGB 
than the AGB (cf. iMcDonald et al.ll2009| ). 

Mass loss can arise from acoustic and/or 
electro-magnetic ch romos p heric activity 



(lHartmann fc MacGregorl Il980t lHartmann fe Avrett] 
119841 : Ivan Loon et al.1 l2010f h or from pulsat io nally- 
enhanced, dust-driv en winds (e.g., iBowenl 119881 : 
Ivan Loon et al.1 120081 ). Radial pulsations levitate ma- 
terial to radii where the stellar atmosphere is cool, 
and pulsational shocks provide the density required for 
dust condensation. Radiation pressure on the grains 
and dust-gas momentum coupling drive the wind. 

1 STScI, 3700 San Martin Drive, Baltimore, MD 21218 USA; 
mboyer@stsci.edu 

2 Astrophysics Group, Lennard-Jones Laboratories, Keele Uni- 
versity, Staffordshire ST5 5BG, UK 

3 Jodrell Bank Centre for Astrophysics, Alan Turing Building, 
University of Manchester, M13 9PL, UK 

4 Department of Astronomy, University of Wisconsin, Madison, 
475 North Charter Street, Madison, WI 53706-1582 USA 

5 Steward Observatory, University of Arizona, 933 North Cherry 
Avenue, Tucson, AZ 85721 USA 

6 Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, MS 65, Cambridge, MA 02138-1516 USA 

7 Department of Astronomy, University of Virginia, P.O. Box 
3818, Charlottesville, VA 22903-0818 USA 

8 Space Science Institute, 4750 Walnut Street, Suite 205, Boul- 
der, CO 80301 USA 



However, it is uncertain whether the dust is the driving 
mecha nism or is simpl y form e d as a byproduct of the 
wind (Iron" Loon et alJ 120051: iFerrarotti fc Gail 120061: 
iHofner fe Andersenll2007t lHofnerll2008h~ 

With the launch of the Spitzer Space Telescope 
([Werner et al.ll2004l ). circumstellar dust has become eas- 
ily detectable, leading to a more complete picture of dust- 
driven mass loss: (1) dust forms only at or above th e tip 
of the RGB (TRGB) (e.g.. IBover etai1[2008l [2009al) : (2) 
dust formation is not inhibited at low me tallicity (e.g., 
IBover et al.ll2006L l2009bHSloan et al.ll2009f ): and (3) dust 
formation a ppears episodic, pos s ibly tied to the pulsatio n 
cycle (e.g.. lOriglia et al.l l200l IMcDonald et al.l l2009f) . 
The first of t hese points has been more controversial. 
iKalirai et al.l ([2007) speculate that strong mass loss in 
RGB stars might explai n the presence of He wh ite dwarfs 
in NGC 6791, although Ivan Loon et all (|2008l ) find that, 
along with an absence of dust, the RGB and HB lumi- 
nosity function is not depleted. lOriglia et al.l (|2007D find 
93 stars with infrared (IR) exce ss as faint as the HB 
in 47Tuc, and IBover et al.l {2006) find 12 such sources in 
M15 (see their Fig. 3b). The presence of dusty stars so far 
down the RGB is surprising given low luminosities, warm 
photospheres, and a lack of pulsation. Indeed, more re- 
cent Spitze r observations show very few d usty stars in 
NGC 362 (IBover et al I l2009aT) and u Cen (|Bover et al.l 
l2008t IMcDonald et al.ll2(J09TV 

This letter examines dust production in 47 Tuc, which 
is am ong the most massive (My = —9.42 mag 9 ; [ Harris 



fl99l ). metal-rich ([Fe/H] = -0.7: lThompson et al II2009F 
and nearby GCs ((m-M) = 13.32; lFerraro et al.lll999fl . 
We present ne w Spitzer observatio ns along with a re- 
analysis of the lOriglia et al.l (|2007f) data. We discuss 
how stellar blending affects the low-resolution 8 /im pho- 
tometry and find that there is little, if any, dust below 
the TRGB in 47 Tuc. 



2. DATA AND ANALYSIS 



9 The Harris (1996) catalog was updated in 
http: / / www.physics.mcmaster.ca/~harris/mwgc.dat 
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Fig. 1. — DSS image of 47 Tuc showing the Spitzer coverage. The 
SAGE-SMC coverage used here is shown in black and the Rood 
coverage in white. The 3.6 and 8 fira coverages are represented by 
dashed and solid lines, respectively. We limited the SAGE-SMC 
coverage to R < 10' to minimize contamination from the SMC and 
to mimic the Rood coverage. 

The small separation («2°) between 47 Tuc and the 
Small Magellanic Cloud (SMC) resulted in serendipitous 
Spitzer observations of the south-eastern half of 47 Tuc 
(Fig. [I]) as part of the Spitzer Surveying the Agents of 
Galaxy Evolution Legacy program (SAGE-SMC). The 
SAGE-SMC observations consist of two epochs, sepa- 
rated by three months, of Infrared Array Camera (IRAC) 
maps at 3.6, 4.5, 5.8, and 8 /im. Here, we use photometry 
obtained by combining both epochs. Exposure times are 
short enough (12 s) to ensure that bright sources remain 
unsaturated, but long enough to obtain sufficient sen- 
sitivity for reasonably comple te photometry to beyond 
the HB (cf. iBover et al.ll2009a[) . Angular resolution with 
IRAC ranges from 1.7" at 3.6 /im to 1.9" at 8 /im. The 
SAGE-SMC MIPS observations (24 - 160 (jm) did not 
cover 47 Tuc. For details regarding data acquisition and 
reduction for SAGE-SMC, see K.D.Gordon et al. (2010, 
in preparation). 

47 Tuc was also observed with IRAC as part of 
PID20298 (P.I. R.Rood). We obtained these data from 
the Spitzer archive and reduced them independently 
(hereafter referred to as the Rood data). The Rood 
field-of-view at each IRAC wavelength reaches a clus- 
ter radius of «10', with an overlapping area of 5' x 9' 
centered on 47 Tuc (Fig. [IJ. Two depths were achieved: 
12 s pix -1 and 936 s pix -1 , hereafter described as the 
shallow and deep da ta, respective l y. O bservation de- 
tails are described in lOriglia et al.l ([2007f ) . Here, these 
data were processed with pipeline S18.7.0. The Basic 
Calibrated Data were reduc ed and mosaicked with the 
MOPEX reduction package (jMakovoz fc Marleaul [20051) 
after applying an array distortion correction. We im- 
plemented the MOPEX overlap routine to match back- 
grounds between overlapping areas of the images and the 
MOPEX mosaicker for outlier elimination, image inter- 
polation, and co-addition. The final mosaics have pixel 
sizes of 0.86" pix -1 . 



3. PHOTOMETRY 

K.D.Gordon et al. (2010, in preparation) contains de- 
tails about SAGE-SMC point-source photometry. 

Point-source extraction (> 4<r) was pe rformed on th e 
Rood data with the DAOphot II package (Stetson 1987), 
with point-spread functions (PSFs) derived from >10 
isolated stars in each mosaic. Extended sources and 
outliers broader or narrower than the PSF were re- 
jected based on sharpness/roundness cut-offs. The fi- 
nal fluxes are color-corrected using a 5000 K blackbody 
\Svitzer Science Center! [20061) . an appropriate tempera- 
ture for a typical RGB star. This correction differs by 
<1% from that for a 2000 K blackbody. A pixel-phase- 
dependent correction (|Reach et al.l 120051 ) was applied to 
the 3.6 /im photometry. Photometric errors include stan- 
dard DAOp hot errors and th e IRAC absolute calibra- 
tion errors (jReach et al.l 120051 ). Magnitudes relative to 
a Lyr (Vega) are derived using th e zero magnitudes from 
the S pitzer IRAC data handbook \Svitzer Science Center! 
2006). False star tests including both the core and out- 
skirts of the cluster indicate the photometry is >85% 
complete at the HB in both the deep and shallow mosaics 
at all four IRAC wavelengths. Bolometric corrections 
were det ermined using 2MAS S J- and i\~ s -band pho- 
tometry flS krutskic et al. 2006|), the transformation from 
IMontegriffo et al.l ([1998( 1 and the reddening (E(B-V) = 
0.04 m ag) and distance mod ulus ((m— M)q — 13.32 mag) 
used in lOriglia et~aTl ([2007D . 

The shallow Rood data are confusion limited at R < 2' 
(Fig. [2]). Our analysis suggests that the 3.6 /im data 
deviate from a K ing profile at w2.2' (r corc = 24"; 
IHowell et aL |[2000|), while the 8 / im da ta deviate near 
1.8'. Assuming the lOriglia et al.l ([20071 ) analysis is lim- 
ited by the 8 /im data (since they use high-quality K- 
band data instead of 3.6 /im), we might expect to have 
found up to 25% fewer sources in the inner 2' than were 
detected in their analysis due to greater incompleteness 
at 3.6 /im compared to 8 /im. 

The SAGE-SMC 3.6 and 8 /im mosaics cover the south- 
eastern side of the cluster with 3' < R < 10'. We in- 
clude only sources within 10' to match the Rood coverage 
(Fig. [IJ and limit contamination from the SMC, which 
begins to dominate the source density at R >15' (Fig. [2]). 

4. DISCUSSION 

4.1. Dust Excess 

The CMDs in Figure [3] include SAGE-SMC and Rood 
data. RGB star bolometric magnitudes peak in the near- 
IR, and warm circumstellar dust is detectable as photo- 
metric excess at 8 /im. The envelopes we attempt to 
detect are optically thin at IR wavelengths (those that 
are not would stand out); attenuation of starlight can 
thus be neglected, whilst dust emission is insignificant at 
3.6 /im. Hence it does not mat t er wh ether the K — [8] 
color is used, as in lOriglia et al.l (j2007t ). or the [3.6] — [8] 
color. Indeed, in the re-reduced Rood data, all stars with 
red K s — [8] colors on the upper half of the RGB (where 
photometric scatter is small) also have red [3.6] — [8] col- 
ors. We also note that the mean K s — [3.6] sage color 
for 47 Tuc is a mere 0.08 mag, insignificant compared to 
the range in dust e xcess (0 < K — [8] < 1) suggested by 
lOriglia et all (|2007l) . We prefer to use the [3.6] - [8] color 
because it does not mix different origins of data, and we 



•RGB Dust Production in 47 Tuc' 



3 




0.1 1.0 10.0 

Cluster Radius (arcmin) 



Fig. 2.— Top: Source density profile in SAGE-SMC data. The 
SMC begins to dominate beyond R > 15'. We restricted SAGE 
photometry to R < 10' (see text). Bottom: Source density profiles 
for the shallow Rood data. Dashed lines are King profiles with 
Rcovc = 24". Note the data deviate from the King profiles due to 
source confusion, inside 2.2' at 3.6 /im and 1.8' at 8 fim. 

can confidently apply the stellar blending test described 
in Section l4~4l 

4.2. SAGE-SMC CMD 

The SAGE-SMC CMD is presented i n Figure [3 
The b lue line represents an isochrone from iMarigo et al.1 
(2008) with the appropriate age (11.3 Gyr) and metallic - 
ity (Z = 0.0024) for 47 Tuc from lThompson et all £2009); 
the solid and dashed lines represent the RGB and AGB, 
respectively. The isochrone indicates that the TRGB is 
located at Mboi ~ —3.5 mag, agreeing with the lumi- 
nosity markin g the onset of dust pr oduction in 47 Tuc 
(L w 2000 £ ? - iLebzelter et alJl200(l . 

Aside fro m 3 variable stars wi thin 0.5 mag of 
the TRGB (|Lebzelter fc Wood! 12001 . only one source 
wl mag brighter than the HB is red, with [3.6] — 
[8] w 0.25 mag (R.A.= 00 h 25 m 39M, decl.= -72°08'21'/0, 
J2000). This source is not identified as variable by e ither 
IWeldrake et all (f200l or lLebzelter fc Woodl (p005h . but 
its 3.6 /im magnitude varies by 0.14 mag between the 
two SAGE epochs. This source is isolated, and is not 
affected by stellar blending or other imaging artifacts, so 
the IR excess is likely real. However, the source is in the 
outskirts of the cluster, in the direction of the SMC, and 
may not belong to 47 Tuc. If at the distance of the SMC, 
its absolute 3.6 /im magnitude is M^.q — —7.75 mag, 
which is > 1 mag brighter than the TRGB and con - 
sistent with a typical AGB star (jBolatto et all 120071) . 
Thus we conclude that significant IR excess in 47 Tuc 
is limited to < 1 mag from the TRGB in the SAGE- 
SMC data, simi lar to what is seen in qjCen, NGC 6791 , 
and NGC 362 dBover et~all[200l Ivan Loon et aH 120081 : 
iBover et al.ll2009aD . 

4.3. Spitzer Archive CMDs 



Photometry extracted from the Rood data is notice- 
ably noisier than the SAGE-SMC photometry (Fig. |3l 
middle, right panels). This could partially be due to dif- 
ferences in data acquisition and reduction, but it appears 
the scatter is caused mainly by inclusion of the inner 3' 
of the cluster (light gray dots in Fig. |3|), where crowding 
is an issue. The deeper photometry suffers from satura- 
tion at 3.6 /im already more than a magnitude below the 
TRGB. For clarity, saturated sources are excluded from 
the CMD in Figured 

We see immediately that the new photometry has 
resulted far fe wer I R-excessive sources than found by 
lOriglia et all ([2007I ). That study found > 90 IR- 
excessive stars (R < 2') distributed evenly along the 
entire RGB, while we see only a handful of red sources 
only near the upper RGB, or in increasing number and 
redness towards magnitudes fainter than «2 mag below 
the TRGB. In the shallow data {R < 10'), we find 24 
sources with > 3a IR excess ([3.6] — [8] > 0.17 mag) 
and Mboi < ma g, four times fewer than found by 
lOriglia et ail (|2007t) in the same area of their CMD. This 
is consistent with a n AKARI analysis of 47 Tuc presented 
bv llta et al.l (|2007| ). where it can be seen in their Fig. 2 
that most dusty stars reside near the TRGB. 

Thirte en of these 24 red s t ars ar e variable stars iden- 
tified bv ILebzelter fc Woodl (|2005D . The brightest vari- 
ables, V26 and LW10, may be AGB stars, whilst the 
others trace the RG B. Curiously, both sho w no excess in 
8 — 13 /tm spectra (|van Loon et al.l 12006). but do show 
moderate excess in the Rood data ([3.6] — [8] « 0.2 mag), 
which were taken one month prior to the spectral obser- 
vations. The difference could be due to episodic dust 
production during the pulsation cycle. The variables 
VI and V8 are the reddest stars in the shallow Rood 
data, ag reeing with the dust excess seen in their mid-IR 
spectra (iRamdani fc Jorissenl [20011: lOriglia et al.1 l200l 
ILebzelter et al.ll2006l : Ivan Loon et al.l 120061) . 

When excluding dusty variable stars, we find only 11 
sources with IR excess in the shallow Rood data, all but 
one more than w2 mag fainter than the TRGB. These 
sources all fall within 3' of the cluster center, suggesting 
that stellar crowding may be leading to inflated fluxes, 
as measured in the lower resolution 8 /im images. 

In the deep data, we find 11 other sources with [3.6] — [8] 
excess > 3 a. Ho wever, exclud i ng th e inner 2' of the 
cluster as done in lOriglia et al.1 (|2007l ) leaves only one 
source with IR excess. 

4.4. Stellar Blending and Imaging Artifacts at 8 /Jtm 

A re-analysis of the photometry has revealed only 24 
red sources, missin g 74% of the 93 red sources from 
lOriglia et al.l (|2007| ). This discrepancy cannot be en- 
tirely due to the 25% incompleteness of the 3.6 /tm 
data described in Section [3] since that would require 
red sources to be preferentially undetected compared to 
non-red sources. That scenario is unlikely since the red 
sources are evenly distributed in the cluster and since 
optically thin circumstellar dust emission is minimal at 
3.6 /im. To investigate whether the 24 red sources are 
truly harboring circumstellar dust, we have performed 
a stellar blend i ng te st, similar to that done for to Cen 
in lBover~et al. (2008|). We convolved the 3.6 /tm image 
with the 8 /im PSF and extracted point-source fluxes 
from the resulting low-resolution 3.6 /im image using the 



4 



Boyer et al. 



SAGE-SMC 



Rood Data (shallow) 



Rood Data (deep) 



i i I i i i i I i i i i | i i i i | i i i i | i i 



i I i i i i I i i i i I i i i i | i i i 



T 



i i i i i i i i i i i i i i i i 



-4 " 



-2 " 



- 



2 " 




V23 



SMC AGB star? 
; (see text) 



LW10. 



V26 




' ' ' ' ' ' I ' ' ' ' ' ' ■ ' ' ' ' ' ' 




o o 



R < 3' ■ 
3' < R < 10' o 

Variable Stars: 
R < 3' + 
3' < R < 10' * 



-0.5 0.0 0.5 1.0 1.5 

([3.6] - [8]) 



-0.5 0.0 0.5 1.0 
([3.6] - [8]) 



1.5 -0.5 0.0 0.5 1.0 

([3 6] - [8]) 



1.5 



Fig. 3. — CMDs of47Tuc. Lef t : SAG E-SMC, Middle: Shallow Rood data, Right: Deep Rood data. Note that we see far few er red sources 
here than seen by Origlia et al. (2007) also in the same original data. The blue line in the left panel is an isochrone from lMarigo et alJ 



pOOl) , with Z = 0.0024 and t a 



11.3 Gyr; the dashed portion marks the AGB. The black dot-dashed line marks the TRGB luminosity, 



as indicated by the isochrone. Saturated sources are excluded from the deep CMD. 
while the other variable stars trace the RGB. 



Variable stars LW10 and V26 may trace the AGB, 



same procedure described in Section [3j In the ideal sce- 
nario, where all stars are isolated, [3.6] rca i — [3.6] convo ivod 
should always be near zero. Any sources that show 
[3. 6] rca i — [3.6] convolved > mag are affected by stellar 
blending, which inflates the convolved 3.6 /im flux. Stars 
that are truly red should have [3.6] rca i — [8] > mag 
and [3.6]rcai — [3. 6] convolved ~ mag. Figure 0] shows the 
results of this test. 

In the crowded environment of a GC, blending can be 
complicated. The test described above may also result in 
blue colors, caused by stellar confusion within compact 
groups of stars. With many overlapping PSFs, it is no 
surprise that flux measurements can be uncertain. We 
choose a conservative cut-off of | [3.6] rca i — [3. 6] convolved > 
0.05 mag to define stars with affected photometry. 

The boxes in Figure 0] show the mean 3<r [3.6] — [8] ex- 
cess in the x-direction and the cut-off for blended sources 
in the y-direction. The shaded regions to the left and 
right of these boxes show the approximate locus of truly 
red or blue IR colors. Stars above or below the shaded re- 
gion are suffering from blending and/or source confusion. 
Stars shown to have real red or blue [3.6] — [8] colors are 
marked in red and blue, respectively, and sources affected 
by blending and/or confusion are marked in green. The 
locations of these stars on the CMDs are shown in the 
lower panels. In the shallow data, the horizontal plume 
of stars at the TRGB is striking and very well distin- 
guished from the vertical dispersion, which is a result of 
stellar blending and source confusion. Note that in the 



shallow data, only one red star with Mboi ^ —3 mag is 
blended: variable LW8. 

Of the II sources in the shallow data (identified in 
Section |4~B1 brighter than the HB and not identified as 
variable stars, we find a total of 6 unblended, red sources 
below the TRGB, all with M bo i > -2 mag. In the deep 
data, we find 5 such sources. Tellingly, none of these 
sources arc red in both the shallow and deep data. 

Closer inspection of the shallow 8 /im mosaic shows 
that all 6 of the faint, red, unblended stars fall di- 
rectly on image artifacts caused by internal optical scat- 
tering inside the detector array (known as banding; 
Spitzer Science Center! [2006h . causing an artificial in- 
crease in their 8 /im fluxes. Therefore, there are no truly 
red sources below the TRGB in the shallow data. Of 
the 5 faint, red, unblended stars in the deep data, 2 are 
affected by banding, and 2 lie immediately adjacent to a 
saturated star, which may affect both their 3.6 and 8 /an 
fluxes. In the deep data, we therefore find only one po- 
tentially dusty source below the TRGB. Again, we note 
that this source is not red in the shallow data, indicat- 
ing that the IR excess is not real. Therefore, there are 
no truly dusty sources fainter than wO.5 mag from the 
TRGB, consistent with the SAGE-SMC data. 



4.5. Implications for Mass Loss and Dust Production 

To b e consistent with the analysis from lOriglia et all 
(|200l . we would have detected up to 70 dusty sources, 
even after considering incompleteness in the 3.6 /im data 
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Fig. 4. — Upper panels: Sour ce blending and confusion at 8 /im 
above the HB. See Section I4.4I for a description of the test. Lower 
panels: CMDs showing the locations of stellar blends, truly red 
sources, and truly blue sources. The mean 3<r excess ([3.6] — [8] > 
0.17 mag) is marked by a dotted line. Note most red stars are 
confined to within 1 mag of the TRGB, and most stellar blending 
occurs among the faintest stars. 

compared to their K-band data. However, through close 
inspection of stellar blending and imaging artifacts in 
SAGE-SMC data and the re-reduced lOriglia et al.l (|20Q7f ) 
data, we find that the number of truly dusty stars in 
47 Tuc below about a magnitude fainter than the TRGB 
is consistent with zero. We therefore find that the mass 



loss law derived by IQriglia et al.l J2007) over-predicts 
mass loss along the RGB - estimating nearly two orders 
of magnitude more mass loss among the faintest RGB 
stars than Reimers' law ([Reimersl [1975) and attribut- 
ing «25% of the total mass loss to stars fainter than 
Mho\ = —1-5 mag. With no dusty mass loss along the 
RGB, these faint RGB stars together contribute closer to 
1% of the total cluster mass loss (based on Reimers' law). 
We do see a large number of red sources at or above the 
TRGB (Fig. [3]), as expected if the stellar structure only 
becomes favorable to dust production if the star is either 
close to the occurrence of core helium ignition (i.e., near 
the TRGB), or is experiencing helium-shell flashes (i.e., 
on the thermally pulsing AGB). 

Very little, if any, dust is forming below the TRGB, 
but that is not to say that no mass is lost on the RGB. 
Indeed, the morphology of the HB requires that mass 
is lost on the RGB. Chromospherically-driven winds 
are indicated in several RGB stars by asymmetries and 
coreshifts in chromospheric spectral lines, with mass- 
loss rates up to 10~ 8 M P) yr" 1 (e.g.. IMauas et all [2001 
IMcDonald fc van Looi]||2007t IDupree et al.H2009f ). which 
is almost exactly what is required to explain the HB mor- 
phology. 

In addition to only residing near the TRGB, the ma- 
jority of the stars showing evidence for circumstellar 
dust are also variable. This supports the idea that dust 
production coincides with pulsationally-enhanced stellar 
winds. 
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